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A sterile neutrino of ~ keV mass is a well motivated dark matter candidate. Its decay generates an 
X-ray line that offers a unique target for X-ray telescopes. For the first time, we use the Gamma-ray 
Burst Monitor (GBM) onboard the Fermi Gamma-Ray Space Telescope to search for sterile neutrino 
decay lines; our analysis covers the energy range 10-25 keV (sterile neutrino mass 20-50 keV), which 
is inaccessible to X-ray and gamma-ray satellites such as Chandra, Suzaku , XMM-Newton, and 
INTEGRAL. The extremely wide field of view of the GBM enables a large fraction of the Milky 
Way dark matter halo to be probed. After implementing careful data cuts, we obtain ~ 53 days of full 
sky observational data. We observe an excess of photons towards the Galactic Center, as expected 
from astrophysical emission. We search for sterile neutrino decay lines in the energy spectrum, and 
find no significant signal. From this, we obtain upper limits on the sterile neutrino mixing angle as 
a function of mass. In the sterile neutrino mass range 25-40 keV, we improve upon previous upper 
limits by approximately an order of magnitude. Better understanding of detector and astrophysical 
backgrounds, as well as detector response, will further improve the sensitivity of a search with the 
GBM. 

PACS numbers: 95.35.+d, 13.35.Hb, 14.60.St, 14.60.Pq 


I. INTRODUCTION 

Right-handed neutral fermions (henceforth sterile neu¬ 
trinos) arise in many extensions of the Standard Model in 
explaining the observed flavor oscillations of active neu¬ 
trinos, and yield an extremely rich phenomenology (for 
recent reviews, see, e.g., mi Sterile neutrinos may be 
produced in core-collapse supernovae [5], providing a new 
mechanism for explosion [5], and may explain the origin 
of strong neutron star kicks EMj. The sterile neutrino 
can modify big bang nucleosynthesis CHS] , assist reion¬ 
ization CMS, and affect neutrino oscillations [16j . 

Moreover, it has been noted that sterile neutrinos can 
contribute the entirety of the observed dark matter den¬ 
sity. They could be produced in the early universe via 
oscillation mechanisms, including non-resonantly EZ! or 
resonantly with active neutrinos [18] , or alternatively via 
non-oscillation mechanisms, such as decays of heavy par¬ 
ticles (see mm\ for some of the scenarios). Sterile 
neutrinos could be a warm or cold dark matter candi- 
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date P5H29) . In addition, some sterile neutrino dark 
matter models can explain the baryon asymmetry in the 
Universe [30H55] . 

For sterile neutrinos produced via oscillations to be a 
viable dark matter candidate, they typically have mass 
in the 1 - 100 keV range [TT] UHl US]. They can radia- 
tively decay into an active neutrino and a photon [331 135] . 
The photon carries half of the total energy, and there¬ 
fore lies in the X-ray energy range. The photon line is 
strongly distinct from most astrophysical and detector 
backgrounds, which have smooth energy spectra. An ex¬ 
ception is line emissions from hot gases and activated de¬ 
tector materials. While the decay lifetime must be com¬ 
parable to the age of the Universe in order to ensure that 
sterile neutrinos remain a viable dark matter candidate, 
the decay in high concentrations of dark matter, e.g., cen¬ 
ters of galaxies, clusters of galaxies, and dwarf spheroidal 
galaxies, can lead to an appreciable X-ray flux. The large 
expected flux in many targets, coupled with the spec¬ 
tral and morphological characteristics of the signal, make 
searches with X-rays a very powerful approach for testing 
sterile neutrino dark matter scenarios (see [33UM1I3Z] for 
a comprehensive discussion of various searches). 

The X-ray constraint on sterile neutrino dark mat¬ 
ter was first obtained using the Cosmic X-ray Back¬ 
ground (CXB) [38]. Subsequently a more detailed anal¬ 
ysis was carried out in [39], where the authors consid¬ 
ered a set of galaxy clusters, two spiral galaxies, and 
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the Cosmic X-ray Background (CXB), and used observa¬ 
tions by Chandra and XMM-Newton. Since then, a host 
of other sources have been explored (for a full discussion, 
see, e.g., review articles [331ES]), including other galaxy 
clusters such as Coma [3UJ HI] , the distant A520 1 32] , 
and the Bullet [IB]; nearby galaxies such as Andromeda 
[43M7j and M33 |48| : additional analysis of the CXB 
[T5T - fo7] ; more recently the Milky Way satellites including 
the Large Magellanic Cloud m, Ursa Minor E2 152, 
Draco [55], Willman I [56], and Segue i EH; and, finally, 
the nearest dark matter concentration, the Milky Way 

galaxy [50115111551I55H5T] . 

Several possible line detections consistent with ster¬ 
ile neutrino dark matter decay have been reported. A 
line feature was observed in Willman I that could be 
interpreted as the decay of a sterile neutrino of mass 
m s ~ 5 keV and sin 2 0 « 1CU 9 , where 0 is the mixing 
angle between the sterile and active neutrinos [5'6j (but, 
see [BO’, 63] )■ In another study, X-ray line ratios showed 
an excess that could be interpreted as arising from the 
decays of 17keV sterile neutrinos with sin 2 # ps 10 -12 
m- Most recently, an anomolous X-ray line was de¬ 
tected from galaxy clusters and Andromeda [62 [65] (also 
see [66U73] )- which can be interpreted as the decay of 
7keV sterile neutrinos [22- 



fig. 1: The J-factor, (Eq. Q), as a function of half 

opening angle i/> relative to the GC, for four Milky Way dark 
matter halo profiles. 


II. EXPECTED SIGNAL FLUX 


In this work, we use the Gamma-ray Burst Monitor 
(GBM) onboard the Fermi Gamma-Ray Space Telescope 
to search for X-ray lines. Notable advantages of the GBM 
include its all-sky coverage, which allows the entire Milky 
Way dark matter halo to be explored, and large effective 
area, yielding a very high statistics data set. The energy 
range of the GBM extends from 8 keV up to 40 MeV, con¬ 
veniently filling a gap in energy above the range of pre¬ 
viously considered X-ray satellites and below the range 
of INTEGRAL. Therefore, in this work we focus on this 
unexplored photon energy range E~ = 10-25 keV ( m s = 
20-50keV). We consider the Milky Way because of its 
proximity and well-studied dark matter distribution, and 
because the GBM detectors are more sensitive to large 
scale diffuse emission such as from the Galactic halo, due 
to GBM’s large field of view (FOV) and poor angular 
resolution. 


The primary decay channel of sterile neutrinos is into 
three light active neutrinos. The radiative decay into an 
active neutrino and a photon that we are interested in is 
suppressed by a factor of 27a/87r « 1/128 relative to the 
primary decay channel [35, , and has a decay rate [32 [31] 
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where we have assumed a Majorana sterile neutrino (for 
a Dirac sterile neutrino the decay rate is halved). The 
energy luminosity of decay photons arising from a sterile 
neutrino dark matter clump of mass Mqm is given by 
L 1 = £ 7 (Mdm/ m s )r s , where E 7 = m s /2 is the photon 
energy, and equals 


We describe the expected X-ray signal from sterile neu¬ 
trino dark matter decays in Sec. [TT] The GBM instru¬ 
ment and the dark matter signal modeling in the con¬ 
text of the GBM detectors are presented in Sec. |III| T he 
data reduction procedures are described in Sec. |IV[ In 
Sec. [V] we describe the line search analysis and the pro¬ 
cedure used to obtain limits on sterile neutrino decays. 
We summarize in Sec. m Throughout this work, we 
adopt cosmological parameters from Planck EE where 
H 0 = 100 h km/s/Mpc, h = 0.673, D A = 0.685, n M = 
0.315, h(z) = a/Da + Qm(1 + z ) 3 , the dark matter frac¬ 
tion LIdm = 0.265, and p c = 1.05 x 10 ~ 5 h 2 GeV cm -3 . 
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( 2 ) 

for a typical galaxy-size dark matter halo mass. It can be 
immediately appreciated that this is comparable to the 
total luminosity of astrophysical X-rays in the Milky Way 
in the 2—10 keV range, ~ 10 39 ergs -1 [76], or the total 
Milky Way diffuse emission in the same energy range, 
~ 10 38 ergs -1 [77]. 


The photon intensity (number flux per solid angle) of 
sterile neutrino dark matter decay coming from an angle 
if) away from the Galactic Center (GC) consists of both 
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TABLE I: Dark matter profile parameters for widely adopted 
dark matter profiles in the literature. Our canonical profile is 
the NFW profile. 


Profile 

a 

p 

7 

Rs [kpc] 

NFW 

1 

3 

1 

20 

cNFW 

1 

3 

1.15 

23.7 

Cored isothermal (ISO) 

2 

2 

0 

3.5 

Einasto (EIN) 

- 

- 

- 

20 


the Galactic and the extragalactic components, 
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where t s = 1 /T S is the lifetime, p Q = 0.3GeVcm -3 is 
the local dark matter mass density, R & = 8.5 kpc is the 
Sun’s distance to the GC, and dN/dE = 5(E — m s / 2) 
is the dark matter decay spectrum. The first term in 
the bracket is the Galactic component. The so-called 
J-factor, J(ip), is the integral of the dark matter mass 
density p in the Milky Way halo along the line-of-sight, 


= — 1 ~5~ j d£p(ip,£), (4) 

PqR& Jo 

where £ ma x is the outer limit of the dark matter halo. 
We assume the dark matter distribution is spherically 
symmetric about the GC, hence 


p(ip,£) = p(r G c (tM)) = p [\Jr% - 2 f R q cos ip +1 2 'j . 

(5) 

The value of £ ma x differs depending on the adopted 
halo model, but the contribution to J (ip) from beyond 
~ 30 kpc is negligible. We adopt Umax = 250 kpc in this 
work. 

The second term in the bracket of Eq. ([3]) describes the 
isotropic extragalactic component, where E' = E( 1 + z). 
The factor Reg roughly compares the contribution of the 
extragalactic component versus the Galactic component, 
up to the shape of the energy spectrum. 


Reg = 


c SIdmPc 

Ho PqRq 


~ 2 . 


( 6 ) 


Normally, the extragalactic component can be ignored 
as typically the analysis region is chosen to be a small 
patch of the sky where the Galactic component is much 
larger (e.g., the GC, where J > 1). However, in our 
case, the large FOV of the GBM makes the extragalactic 
component non-negligible. 

The dark matter density profile p(r) of the Milky Way 
is not precisely known, in particular at small Galactic 


radius. We consider several fitting functions that cap¬ 
ture the results of numerical simulations of dark matter 
halo profiles, which can be parameterized by the follow¬ 
ing form, 


p aM {r) = p & 


( r V 7 

1 + (Rq / R s ) a 

\Ro) 

1 + (r/R s ) a _ 
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where parameters for commonly used profiles are summa¬ 
rized in Table |l] Another profile favored by recent simu¬ 
lations is the Einasto profile, 


p Ein (r) = pq exp — , ( 8 ) 

with cxe = 0.17 and scale radius R s = 20 kpc. These 
profiles differ mainly at small Galactic radius. The first 
three profiles have constant logarithmic slopes at small 
radii, which are described by the 7 factor. The Einasto 
profile has the same slope as the NFW profile at the scale 
radius, but the slope decreases as the radius decreases. 

In Fig. |TJ we show the J-factor J (ip) for each dark mat¬ 
ter profile as a function of the angle ip viewed away from 
the GC. The differences between profiles are relatively 
small, because the density p appears linearly in the de¬ 
cay flux (as opposed to in the annihilation flux where the 
density appears quadratically). We use the NFW profile 
as our canonical profile in this work. As will be shown 
in Sec. |IIIB[ the impact of varying the profile is minimal 
after taking into account the detector response and the 
FOV. Thus the sterile neutrino constraint obtained using 
GBM is robust against dark matter profile uncertainties. 

A crude estimate of the expected number of photons 
i / 7 per unit time T from Galactic dark matter decay is 



where we use representative values for the effective area 
and solid angle, the J-factor at ip = 60°, Veo, and a 
nominal sterile neutrino mixing angle. It is immediately 
clear that even a small fraction of the total Fermi-GBM 
live time can yield significant number of signal photons. 


III. INSTRUMENT AND SIGNAL MODELING 
A. GBM Instrumentation 

The GBM consists of 14 detectors: 12 Nal detec¬ 
tors, each operating over energies from 8 keV to 1 MeV, 
and 2 BGO detectors, each operating over energies from 
200 keV to 40 MeV. The Nal detectors are located on the 
corners and sides of the spacecraft, with different ori¬ 
entations, and they together provide a nearly complete 
coverage of the occulted sky. At any given time, typi¬ 
cally 3-4 Nal detectors view the Earth within 60 degrees 
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FIG. 2: The effective area for det-7 Nal detector versus the 
detector zenith angle relative to the detector normal, for three 
example energy bins. Points are data from GBM calibration 
files, and the anomalous dips in the effective area come from 
blockages from other satellite components. Solid lines are fits 
to the data neglecting the dips. 


of the detector zenith, i.e., their FOV is occulted by the 
Earth. 

Not all of the Nal detectors are best suited for dark 
matter searches. At first consideration, det-0 and det-6 
would seem to be the best detectors to use since they are 
aligned close to the LAT zenith (~ 20° offset). However, 
we find that significant parts of the FOV of these two de¬ 
tectors are actually blocked by the LAT itself. Also, half 
of the detectors are pointed towards the Sun all the time, 
and X-ray emissions from the Sun contaminate their low 
energy spectrum. Lastly, some detectors are pointed side¬ 
ways, i.e. ~ 90° relative the LAT-zenith, which suffer 
large FOV blockage from the Earth. Ruling out these de¬ 
tectors, only det-7 and det-9 seem to be suitable, which 
are ~ 45° relative the LAT-zenith. Upon inspection, we 
observe an anomalous spectral feature in the low energy 
spectrum of det-9 compared to other detectors. As a re¬ 
sult, we use det-7 as our fiducial detector for analysis. 
As will be shown below, this analysis is systematically 
limited, rather than statistically limited. Using only one 
detector for this analysis also avoids introducing system¬ 
atic uncertainties from combining multiple sets of data 
from different detectors. 

The Nal detectors have a wide FOV, as seen in 
Fig.® which shows the effective area versus the detector 
zenith angle, 6 , for the det-7 detector. We obtain the 
GBM effective area data from detector response matrix 
files (GS-008). Each file contains the effective area as 
a function of energy, for a specific detector zenith and 
azimuthal angle. In Fig. [2j each point denotes the effec¬ 
tive area extracted versus the corresponding zenith angle 



FIG. 3: The effective area versus the incident photon energy 
for three incident angles. The points are obtained from GBM 
calibration files and the lines are linear interpolations of the 
points. 


from the detector response file. Beyond about 40°, we 
observe anomalous dips in certain azimuthal direction at 
all energies, which is presumably caused by blockages 
from satellite components in the FOV of the detector. 
We remove these anomalous dips by requiring adjacent 
bins deviate no more than ~ 10%. After this procedure, 
the angular dependence of the effective area can be well- 
described by cosine functions, as shown by the solid lines. 

In Fig. ® we show the energy dependence of the ef¬ 
fective area for three representative zenith angles. The 
points are obtained from the detector response matrix 
files, and are chosen from a specific azimuthal angle at 
which the detector FOV is not blocked. Using the cosine 
fits described above, the energy dependence is obtained 
by linear interpolation of the model in energy. Thus, we 
obtain an azimuthally symmetric model of the effective 
area (i.e., it depends only on zenith angle and energy). 
Shown in solid lines in Fig.® are the model for the given 
zenith angles. The FOV blockages slightly reduce the de¬ 
tector sensitivity towards a particular azimuthal angle. 
This effect, however, is not expected to introduce spuri¬ 
ous spectral features in energy spectrum, since the block¬ 
ages affect all energies. We thus neglect these blockages 
and use the smooth angular fits to model the expected 
signal in the next section. 

One important feature of the Nal detectors is that they 
are limited in their photon-tracking capabilities, i.e., one 
cannot simply obtain the photon flux as a function of the 
incidence direction for a specific position on the sky. In 
other words, the tradeoff for the large FOV is poor angu¬ 
lar resolution. Earth occultation techniques can be em¬ 
ployed to obtain photon direction for point source stud¬ 
ies [Z81E2I, but this technique has not yet been demon- 
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strated for diffuse emissions. Fortunately, the lack of 
photon tracking is not very problematic for sterile neu¬ 
trino dark matter decay searches due to the large angu¬ 
lar extent of the expected emission. However, this does 
mean that one cannot accurately construct an intensity 
sky map of the GBM data (Eq. i»- As a result, one 
needs to properly model the signal taking into account 
the detector response to match the observable. In this 
case, the instrumental observable is counts rate (number 
of photons per second), as a function of the Nal detector 
pointing direction. 


B. Expected Signal Modeling 


Given the sterile neutrino decay photon intensity 
we compute the expected number of photons, 
u i : j. for energy bin i and detector sky-pointing direction 
j. The expected number of photons per observing time, 
Tj , from a particular detector pointing direction is then 


dvj,j 

~dT~ 


/ dE dSl{8) / dE 

/^mm J 2n J 

T(V>, E) G (e, £) A eS (E,d) 1, 


( 10 ) 


where E™ ax and E™ ln are the boundaries of the energy 
bin i. We integrate over the hemisphere the Nal detec¬ 
tor points at, i.e., over the detector zenith angle 9, and 
attribute all the photons to pixel j. A position on the 
sky with an angle relative to the GC, t/f, is related to 
the detector zenith angle and the pixel that the detector 
points at through 0 —► ip(0 7 j). The pointing direction 
of the detector is therefore defined by ip(0,j). The fac¬ 
tor G(E , E) takes into account the energy resolution of 
the Nal detector, which we model as a Gaussian with 
width given by the pre-launch calibrations (80', [81]. The 
energy resolution is about 10% for our analysis range. 
And lastly, A e g(E,9) is our Nal detector effective area 
model, which is a function of energy and the detector 
zenith angle, as in Figs. ii 

Using the Dirac-delta function for the energy spec¬ 
trum, the expected signal is 


dfij 

~dT~ 


( m s ) 


PqRq 


rEf 
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Reg / dE / dQ A e g G(E, E) 


ej^-E) 

Eh(E) 


where J ( E,j ) is the “convolved J-factor”, A f{E) is a 
normalizing factor, h(E) = + ^m{iti s /(2E)) 3 , 0 

is the Heaviside step function. We have suppressed the 
argument of A e g for simplicity. 


The convolved J-factor is defined as 

fj(iP)A eS (E,0)dn(e) 


J(E,j) = 


■a m 


( 12 ) 


which takes into account the effect of detector response; 
it represents the J-factor defined by detector pointing 
directions. It depends on the detector pointing direc¬ 
tion through The normalization factor N{E) = 

2irA e g(E, 0) captures the energy behavior of the effec¬ 
tive area. The normalization of this factor is unimpor¬ 
tant as it cancels itself when obtaining dark matter decay 
fluxes/limits. In Fig. |4j we compare the convolved J- 
factor with the normal J-factor defined in Eq. (4]). Once 
the detector response is taken into account, the differ¬ 
ence between different profiles decreases drastically even 
for pointing directions very close to the GC. E.g., for 
10-11 keV bin, the difference in the convolved J-factor 
between NFW versus EIN and ISO is < 1%. Therefore, 
systematic uncertainties due to the choice of dark matter 
profile are minimal. 

In the left column of Fig. [5] we show the modeled 
dark matter maps for the Nal detector from the Milky 
Way halo for several energies. We pixelate the sky 
into 768 pixels of equal solid angle using the HEALPix 
scheme 1 , i.e., each pixel corresponds to a solid angle of 
Af 1 ~ 1.6 x 10 -2 sr. We use the Milky Way contribu¬ 


tion from Eq. (11), which takes into account detector 


energy and angular response. The extragalactic compo¬ 
nent only adds a constant value to the signal map. We 
choose the line energy to be at the center of the chosen 
energy bin. The decay rates for the sterile neutrino sce¬ 
narios are chosen to approximately match the count rates 
of the corresponding data maps (right column, described 
below). By construction, the Milky Way dark matter 
contribution is spherically symmetric, and the large an¬ 
gular extent of the signal is due to the large FOV and 
poor angular resolution of the Nal detectors. 


IV. DATA SELECTION AND REDUCTION 

In this section, we describe the data reduction proce¬ 
dures to improve the data quality and the cuts designed 
to reduce various backgrounds. At the end we obtain 
a data set that can be compared to Eq. to obtain 
limits for sterile neutrino dark matter decay. 

We use GBM daily data from 12-AUG-2008 to 31- 
DEC-2012, a total of 1601 days. We use the CSPEC 
data (GS-002) with nominal 4.096 s time resolution and 
128 channels in energy from 5 to 1402 keV (the first and 
last few energy bins are not usable). We then devise 
several cuts to improve the data quality. The goal is to 
obtain a data set that is representative of the diffuse sky 
emission as observed by the GBM Nal detectors, while 


1 http://healpix.jpl.nasa.gov m 
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FIG. 4: The convolved J-factor (Eq. (121) versus the opening 
angle with respect to the GC, defined by where the 

detector normal is pointing. The difference between different 
profiles is drastically reduced. A small energy dependence 
is introduced from the effective area. The vertical dotted 
line denotes the boundary of our ROI. Shown in grey are the 
theoretical J-factors from Fig. [T| 


minimizing various types of backgrounds. The most dom¬ 
inant source of background is due to cosmic rays interact¬ 
ing with the satellite, directly activating the detector or 
triggering the detector through delayed radioactive de¬ 
cays of the satellite material. 

To this end, we employ the following cuts: 

• LAT cut. 

We first select data sets that are suitable 
for analysis using data flags from Fermi-LAT 
weekly photon files: LAT_C0NFIG=1, LAT_M0DE=5, 
DATA_QUAL=1, R0CK_ANGLE<50, SAA=F. The first 
three conditions ensure the detector configuration 
and data quality are suitable for scientific analy¬ 
sis. The fourth condition ensures that the Earth is 
not in front of the LAT’s FOV, which is approxi¬ 
mately, but not exactly, the FOV of the Nal detec¬ 
tor (we address this in the Earth cut below). The 
last condition excludes the times when the satellite 
is inside the South Atlantic Anomaly (SAA), where 
the high cosmic ray activity significantly increases 
the radioactivity of the satellite. The GBM detec¬ 
tors are turned off during SAA passage, hence the 
observed counts are zero in these time periods. 

The LAT cuts alone, however, are insufficient for re¬ 
ducing background events, because of the different phys¬ 
ical locations of the detectors on the satellite, different 
backgrounds, and the different technologies of the LAT 
and the GBM. We therefore develop new cuts specifically 
for the GBM. 


• Transient sources cut. 

This cut removes the epochs when the GBM detec¬ 
tors detect transient sources, such as gamma-ray 
bursts, direct cosmic-ray hits, solar flares, Galactic 
X-ray transients, and magnetospheric events, etc. 
Though these transients only occupy a small frac¬ 
tion of the observation time, some of them can be 
bright enough to cause the data acquisition system 
to overflow. 

• Extended SAA cut. 

The LAT cut does not completely remove events 
due to passages of SAA. This is because the satel¬ 
lite is intensively bombarded by cosmic rays during 
each passage through the SAA, leaving the satel¬ 
lite in a highly radioactive state even after leaving 
SAA. This effect is even more pronounced for con¬ 
secutive passes through the SAA. In this case, there 
is insufficient time for the satellite to return to its 
normal radioactive state. As a result, orbits passing 
through the SAA consecutively induce anomalously 
high photon count rates even when the satellite is 
outside the SAA. We therefore apply cuts to remove 
the data collected between consecutive passages of 
the SAA, in addition to the times that the satellite 
is physically in the SAA, which are eliminated in 
the LAT cut. Removing these orbits is important 
to reduce events originated from cosmic rays. 

• Earth cut. 

Lastly we apply two cuts on the orientation and the 
position of the Nal detector relative to the Earth. 
We first require that the angle between the Nal 
detector normal and the vector directed from the 
Earth center to the satellite to be less than 50°. 
This is to reduce contamination from the Earth 
limb and occultation from the Earth itself. The 
next cut is on the geomagnetic coordinate. The 
high-altitude cosmic ray activity is directly corre¬ 
lated to the Earth’s magnetic field structure. The 
number of observed background events increases 
with geomagnetic latitude. To minimize this con¬ 
tamination, we select data only when the geomag¬ 
netic latitude is less than |20|°. 

In Fig. [6j we show an example of the data and the cuts 
we adopt to improve the data quality. The data points 
are the counts rates on 20th December 2008 observed 
by det-7. Each dot corresponds to count rates measured 
over ~ 4 s. We select the energy range from 344 keV to 
471 keV, where the data is dominated by the cosmic-ray- 
induced background. 

The first feature that can be seen in Fig. [6] is the se¬ 
ries of epochs with no count rate. This is because the 
detector was shut down when the satellite is in the SAA. 
These epochs are removed in the LAT cut. It is also clear 
that the count rates are anomalously high even after the 
satellite leaves the SAA (i.e., right after the gap), due to 
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10-11 keV, Sm~2(3Ttieta): 6.4E-11 


10-11 keV 



0.0 61.0 [ ph/sec ] 

16-17 keV, Sin~2(2Theta): 2.6E-12 



0.0 74.1 [ ph/sec ] 

40-42 keV, Sin-2(2Theta): 1.2E-14 
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FIG. 5: (Left) Simulated dark matter counts rate maps in Galactic coordinates for several line energies, taking into account 
detector response. For each map, the assumed line energy is contained in the energy bin shown and sin 2 2d (as labeled) are 
chosen to approximately match the observed counts in the same bin. (Right) The final counts rate sample from 4 years of 
data from the Nal detector, which corresponds to 4.6 million seconds (~ 53 days) of live time after data cuts. All of the sky 
maps are pixelated into 768 HEALPix pixels. The pixel position corresponds to the pointing direction of the detector normal. 
The grey pixels are where no observing time is registered after the selection cuts. 


the increased radioactivity of the satellite. These epochs 
are removed in the Extended SAA cut. These two SAA 
related cuts are represented by the red hashed regions. 

The second feature is the the oscillatory shape dur¬ 
ing the middle of the day. Overlaying the data points 
we also plot the location of the satellite in geomagnetic 
latitude (blue line). One can see the count rates are 
correlated with the geomagnetic latitude. We therefore 
remove all the data recorded when the geomagnetic lati¬ 
tude is larger than 1201 °. This cut is represented by the 


grey shaded region, and the removed data points are la¬ 
belled in blue. The choice of a uniform |20|° geomagnetic 
latitude cut is a balancing act between maximizing sky 
coverage and reducing background. More sophisticated 
cuts may be possible. 

The importance of our cuts for improving the data 
quality can be estimated in Fig. [6] The increased count 
rate right after SAA can be a factor of a few higher. Even 
the variation due to geomagnetic latitude can be up to 
a factor of two. Our Extended SAA cut and Earth cut 









Time [ s ] 


FIG. 6: A sample of the raw photon count rates for the GBM-Nal detector from 20-12-2008, from 344 keV to 471 keV, where 
the counts rate is dominated by the cosmic ray induced background. The blue line indicates the location of the satellite 
in geomagnetic latitude (in absolute value). The epochs with no data are when the satellite is physically in the SAA. The 
red hashed regions represent the cuts on orbits that pass through the SAA. The grey shaded regions illustrate the cuts on 
geomagnetic latitude; the data removed by this cut is colored blue. In the end, only black data points in the white regions are 
used for analysis. 


is therefore necessary to reveal the astrophysical compo¬ 
nent, which is comparable to the detector background at 
low energies (shown below). Transient sources does not 
contribute significantly to the total counts, but they can 
dominate a particular sky pixel. Since they only con¬ 
tribute a small fraction of the live time, the Transient 
sources cut is very efficient. 

The final data products obtained are observed counts 
and exposure time over 128 energy bins and 768 sky pix¬ 
els. The total live time of the data product is ~ 4.6 x 10 6 
seconds (~ 53 days). Despite having a huge reduction 
from the raw data, we are still far from statistically lim¬ 
ited, as will be shown below. 

In the right panel of Fig. [5j we show the counts rate 
sky map for the labeled energy bins. At low energies, 
we observe a clear excess towards the GC. We inter¬ 
pret the excess as astrophysical (i.e., non-instrumental- 
related) emissions from the Milky Way. The astrophysi¬ 
cal flux is about ~ 10 -1 cm -2 s -1 if one extrapolate from 
the high energy observations [831185] . which matches the 
observed counts rate of about ~ 10s _1 . The observed 
excess towards the GC also shows a small north-south 
asymmetry, which probably reflects the underlying dis¬ 
tribution of diffuse and discrete X-ray sources. 

For the maps at high energies, the morphology is sig¬ 
nificantly more isotropic than at low energies, with small 
variations that trace orbital structure, as expected from 
cosmic-ray-induced backgrounds. For example, the two 
dark spots near the orbital pole in high energies is also 
seen in the low energy map. 

As a result, we conclude the low energy data set con¬ 
sists of a mixture of astrophysical diffuse and point source 
emissions, plus residual cosmic-ray-induced background. 
The grey pixels in the maps represent positions on the sky 
that were not visited by the detector, and are excluded 


from the analysis. 

Using the data sky map and the convolved J-factor 
we can determine the region of interest (ROI) 
for our analysis. As J(E,j) flattens out at small angles 
due to the poor detector angular resolution, as shown in 
Fig.0 there is little benefit in choosing a small ROI. We 
carry out a signal-to-noise study to look for an optimal 
ROI angle. The morphology of the GC excess seen in 
low energies turns out to be comparable to the smoothed 
dark matter distribution, and the signal-to-noise is fairly 
insensitive to the choice of angle. This is a direct conse¬ 
quence of the poor angular resolution of the Nal detector. 
We conservatively choose a large ROI, which consists of 
pixels within 60° from the GC, i.e., ^ < 60°. With this 
selection, we have enough pixels to average out poten¬ 
tially spurious behavior in some individual pixels, and 
have more than enough statistics. Lastly, this ROI only 
minimally overlaps with the dark spot positions near the 
orbital poles. 

In Fig. [7] we show the binned counts spectrum for the 
data sample in the GC ROI. As a comparison, we also 
show the spectrum for the anti-GC ROI (iJj > 120°). The 
total observed time for the two samples are 975066 s and 
911451 s, respectively, and this difference is the main rea¬ 
son the normalization differ in high energies. In general, 
the counts spectrum has a power-law behavior at high en¬ 
ergies, as expected from cosmic-ray induced background. 
There are several prominent line features from excited 
energy levels of 127 I at 57.6 keV and 202.9 keV, as well 
as the 511 keV line from positron annihilation from the 
atmosphere and nearby materials m- At low energies, 
the GC and anti-GC spectral shape starts to deviate, 
and the difference in normalization increases compared 
to high energies. This indicates that the astrophysical 
component starts to appear in the GC sample. 
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V. LIMITS ON STERILE NEUTRINOS 

We present two limits on sterile neutrino decay lines. 
The first is a conservative limit based purely on flux com¬ 
parison. The second uses the fact that the signal is a 
photon line, while the background flux is approximately 
a power-law within the search energy window. 


A. Flux Analysis 


The most robust constraint one can place on the am¬ 
plitude of a sterile neutrino dark matter decay signal is to 
require that the expected signal counts do not exceed the 
total measured counts. For a set of dark matter masses, 
we compare, bin by bin, the predicted counts from sterile 
neutrino decay to the total counts measured. This ap¬ 
proach therefore assumes the hypothesized signal dom¬ 
inates the observed spectrum without any assumptions 
about the detector and astrophysical background. 

The expected signal counts are given by summing the 
count rates in all the individual sky pixels within the 
ROI, using Eq. weighted by the actual observing 

time Tj in each pixel, 


J'i 


ROI 


E T i 


disij 

~dT~- 


(13) 


The measured photon counts data from all pixels in the 
ROI is 

ROI 

N 'd > ( 14 ) 

i 


where N it j is the number of counts in energy bin i and 
pixel j measured by the GBM detector. 

We obtain the flux analysis limit on the decay rate, 
T s , by requiring i/j < di for all energy bins for each m s . 
The limit obtained this way is very conservative. It is 
unlikely that sterile neutrino decay, which has a sharp 
spectral shape, would dominate a narrow energy range 
in the count spectrum while other components conspire 
to vanish in that particular energy range. 


B. Spectral Analysis 

The sensitivity to sterile neutrino dark matter decay 
can be improved dramatically using the observation that 
the sterile neutrino decay signal and the dominant back¬ 
ground have different spectral shapes. 


1. A simple background model 

As shown above, our GC data sample contains an as¬ 
trophysical component as well as internal detector back¬ 
grounds. The astrophysical contribution from the inner 
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FIG. 7: The counts spectrum for the final data sample for 
both GC ROI {ip < 60°) and anti-GC ROI {ip > 120°) chosen 
to have the same solid angle. The dominant component is a 
power-law plus various background lines. The overall differ¬ 
ence in normalization is due to the higher exposure towards 
the GC than the anti-GC direction. The additional excess at 
low energies towards the GC region suggests the rise of the 
astrophysical component. The vertical dotted line indicates 
the energy bins used for the spectral analysis. 


galaxy is dominated by points sources (all unresolved by 
GBM), and the energy spectrum was shown to be well 
described by a power law above 20 keV [55H85] . The in¬ 
ternal detector background is a consequence of cosmic 
rays interacting with the satellite components, which re¬ 
tains the power-law behavior of the incoming cosmic rays. 
We therefore expect the energy spectrum to have a power 
law distribution. 

A power-law spectrum is an even better approximation 
when we analyze the data in small energy windows. We 
consider 15 of such search windows, one for each line en¬ 
ergy of interest. The line energies are the corresponding 
energies of the energy bin number 6 to 20 in GBM num¬ 
bering scheme (labeled by io). For each search bin i 0 , the 
search window contains a number of energy bins (labeled 
by z), where 

Max (imin, io — A to) < * < zo + Aw. (15) 

The window size is Aw = 5, which makes the window 
width on each side about 3-4 a of the energy resolution 
at the line energy. For line energies near the low energy 
cutoff, we truncate the search window at the lowest us¬ 
able energy bin, i mm = 6, which corresponds to a central 
bin energy of 9.3 keV. The signal line energy in such case 
is not located in the center bin of the search window. 

With the power-law assumption for the non-dark mat¬ 
ter components in each search window, the model pho- 
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ton counts spectrum therefore contains the dark matter 
signal component ( dv/dE ) and a power-law background 
component ( db/dE ), 


dv 

~dE 


db 

dE 


f s ^S(E - E 0 )Af(E) 

ETj ( 


(16) 


-Reg 

P 


ZTjJ(Eo,j) 

E '-' 1 


A cS {E)dn 


0(E o - E) 
Eh(E) 


En 


A r(E), 


(17) 


respectively, where E 0 is the energy of bin i 0 . The factor 
Af(E) takes into account the energy response of the ef¬ 
fective area. The model has only three free parameters, 
fsiPi and 7 . The factor f s is the amplitude of the dark 
matter signal, which is the only parameter that we are 
interested in. The normalization and the spectral index 
of the background power law are thus treated as nuisance 
parameters, 3 = 

The total expected counts in energy bin i in the search 
window is then obtained by convolving with the detec¬ 
tor energy resolution and integrating the model over the 
energy bin, 


Vi + bi= I dE I dE ( ^ + H) G(E,E ). (18) 


Comparing the data model (Eq. (13)) with the expected 
signal (Eq. 01), the line amplitude f s is related to sterile 
neutrino parameters by 


fs = 


PqRq 

Anm s T 




(19) 


„ „ „„_ 2 _2 -1 -1 / sin 2 2 6 

= 8.6 x 10 "cm 2 s 1 sr 1 ( - „_ 1Q 
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J2 T A E o,j)- 


To search for a line signal from the data, it is impor¬ 
tant to understand the uncertainties associated with the 
measurement. We first consider the systematic uncer¬ 
tainty in the effective area of the detector, which is ~ 5% 
according to the GBM collaboration PI 151] . Note the 
quoted uncertainty is the total uncertainty for the effec¬ 
tive area, which in principle can be two different kinds of 
uncertainty. The first kind is the overall uncertainty on 
the effective area across all energy bins, which affects the 
value of the flux obtained from data. The second kind is 
the uncorrelated errors between energy bins, which may 
introduce spurious spectral features even if the true flux 
spectrum and the true effective area are both smooth in 
energy. For a spectral analysis, the uncorrelated error 
among energy bins is much more important than simply 
a normalization shift. In this work, we conservatively at¬ 
tribute all the 5% uncertainty to the uncorrelated errors. 
As a result, the model uncertainty for each energy bin is 


( 20 ) 


We then consider the statistical uncertainty. As shown in 
Fig -0 the number of photons is enormous in the energy 
range that we are interested in. The statistical uncer¬ 
tainty in each bin is small, y/N/N < 10~ 3 . Therefore, 
we safely ignore the statistical uncertainties in this work. 

We adopt the method of maximum likelihood for fitting 
the counts spectrum for each search window. For each 
search bin io, we assume a Gaussian probability distribu¬ 
tion function for each energy bin in the search window, 
giving the likelihood function: 

where the product is taken over the energy bins i in the 
search window. Best fit parameters are obtained by max¬ 
imizing the likelihood function, or equivalently minimiz¬ 
ing its negative logarithm. 

We first find the best-fit background only parameters, 
So, where f s is set to zero. The implicitly defined S 0 is 
given by 

A (f s =0,S o |*o) = Min{- 2 Log£(/ s = 0,E|i 0 );3}. 

( 22 ) 

We check whether the power-law only background model 
is a reasonable hypothesis by computing the reduced 
X 2 (x 2 P er degree of freedom) for each search window. 
We find that the reduced \ 2 ranges from 0.2 to 1.4 in our 
analysis range. We therefore conclude that the power-law 
only model plus the prescribed 5% systematic error can 
reasonably describe the data for each search window. 

In Fig. [ 8 j we show explicitly the 15 search windows 
for this analysis. The blue data points are the GBM 
data, and the assigned error bars are the 5% systematic 
uncertainty. The statistical errors are too small to be 
shown. The blue lines are the best fit count spectrum 
from the power-law only model described above. The 
apparent peculiar spectral features, such as those around 
18 and 26keV, are successfully captured by the power- 
law model when effective area and non-uniform energy 
bins are taken into account. 


2. Limits on dark matter decay rate 

To search or constrain the line signal, we use the so- 
called profile likelihood method i 86 j. We search for the 
best fit line amplitude by minimizing the negative log- 
likelihood with respect to all the model parameters, 


A(/ s o,H 0 |*o) =Min{-2Log£(/ s ,S|* 0 );/ s ,S}, (23) 

where f s is constraint to be non-negative. We observe 
no significant preference for the presence of the line sig¬ 
nal. We then proceed to find the 95% C.L. one-sided up¬ 
per limits on the dark matter signal amplitude, /® 5 , by 
increasing the amplitude while continuously minimizing 
the log-likelihood function over the nuisance parameters, 


CTAeff = 0.05(z/j + bi) . 
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FIG. 8: The blue points are the measured data with error bars indicating the 5% systematic error. The blue line is the best fit 
model to the data with only the power-law component. The red line shows the best fit model when including the line signal 
with 95% upper limit amplitude. The red arrow indicates the central energy of the line signal. 
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FIG. 9: The black solid line is the 95% C.L. upper limit from 
the spectral analysis. The green (yellow) shaded region shows 
the 68% (95%) intervals from the Monte Carlo simulations. 


until it is 2.71 larger than the best fit log-likelihood, 

Min{—2Log£(/f,S|i 0 ); S} = A(/ s0 , E 0 \i 0 ) + 2.71. (24) 

The best fit model when the line signal is at 95% upper 
limit is shown in Fig. [8] in red lines. The red arrow in¬ 
dicates the energy of the inserted X-ray line. Using /® 5 
and Eq. (191, we then obtain the 95% C.L. upper limit 


of the dark matter decay rate. 

We perform a Monte Carlo study to check the robust¬ 
ness of the limit. For each search window, we generate 
100 mock data sets, according to the best fit power-law 
only parameters and the 5% systematic error with Gaus¬ 
sian probability distribution function. We perform the 
profile likelihood analysis to obtain the 95% upper limits 
for the mock data sets. In Fig. |9j we show the obtained 
upper limit from data and the 68% and 95% coverage of 
the limits from our Monte Carlo simulations. 

Overall, we find that the observed limit is consistent 
with our Monte Carlo realizations at the 95% level. At 
about 28keV sterile neutrino mass, we find the actual 
limit touches the 95% lower bound of the expect limit. 
This is likely due to the data point at about 14keV falls 
below what one would expect from a smooth power-law 
flux spectrum, as shown in the fits in Fig. [8} There are no 
known detector defects at this energy |5Tj, we thus 
consider this as a ~ 2er systematic downward fluctua¬ 
tion in the effective area model. This downward shift 
effectively means the data prefers a negative line, which 
results in the improved limit at this energy. It is also im¬ 
portant to note that the limits are correlated due to the 
largely overlapping data points in adjacent search win¬ 
dows. As a result, limits from line energies close to 14 keV 


FIG. 10: The conservative upper limit from the flux analy¬ 
sis (Black) and the 95% C.L. upper limit from the spectral 
analysis (Blue) for the decay rate of sterile neutrino dark mat¬ 
ter. The hashed regions are excluded by the corresponding 
analyses. 


are all slightly improved. The simulated limits from our 
Monte Carlo realizations are not correlated with adja¬ 
cent line energies, since the mock data sets are generated 
independently for each search window. 

Finally, Fig. pH] shows the limits obtained on the decay 
rate from both the flux analysis and the spectral analy¬ 
sis, with the hashed region corresponding to the excluded 
parameter space. As expected, the spectral analysis pro¬ 
duces a much stronger limit than the flux analysis. Since 
the presence of a line signal mostly only affects one energy 
bin, one would expect the spectral analysis limit is ap¬ 
proximately given by the size of the error bars of the data 
points, and thus the spectral analysis limit is expected to 
be about 5% of the flux analysis limit. This is indeed the 
case in most of the mass range, except where the data 
prefers a negative line, as discussed above. At low ener¬ 
gies, the spectral analysis limit deteriorates rapidly. This 
is due to the imposed lower cutoff of the search energy 
window. As the line energy approaches the boundary, the 
number of bins used for the fit is reduced accordingly, and 
the spectral shape becomes increasingly degenerate with 
the power-law shape. Both factors cause the limit from 
spectral analysis to deteriorate. 


C. Limits on Sterile Neutrino Dark Matter 

Using the upper limits on the decay rate, we derive the 
corresponding upper limit on the mixing angles for sterile 
neutrino dark matter. In Fig. El we show the constraint 
on the mixing angle-mass plane. For comparison, we 
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also show the only limit in this energy range, obtained 
with CXB observations using HEAO-1 [50] . and with 
Milky Way observations using INTEGRAL [BO]. Unsur¬ 
prisingly, the flux analysis does not yield competitive lim¬ 
its. However, in the mass range m s ~ 25 - 50keV, the 
spectral analysis improves the limit on the mixing an¬ 
gles by about an order of magnitude compared to the 
previously strongest limit. 

Compared to the previous analysis m, which uses 
HEAO-1 A4 Low Energy Detector data EZI, our analy¬ 
sis improves mainly in: The GBM data has smaller error 
bars compared to the HEAO-1 data (~ 5% vs ~ 10%); 
We have employed several cuts to reduce cosmic rays in¬ 
duced backgrounds; The Milky Way halo yields a larger 
signal flux than the CXB alone; The GBM Nal detector 
has a slightly better energy resolution. These factors, 
each expected to give a factor of few improvement, all 
contribute to our improved limit. 


VI. DISCUSSION AND CONCLUSION 


A. Future Developments 


In this work, we obtain competitive limits on sterile 
neutrino dark matter decay by analyzing GBM data. 
This is the first time the GBM data is used for a dark 
matter search, and we have obtained the strongest con¬ 
straint available in the mass range 25 — 50keV (Fig. 111. 


Although we focus on the implications for the sterile neu¬ 
trino, our limits can be applied to all dark matter candi¬ 
dates that produce a mono-energetic photon in the keV 
range, such as moduli dark matter, gravitino dark mat¬ 
ter, and other candidates [88H90] . It is straightforward 
to constrain the parameter space for the corresponding 
dark matter candidates using the limit on the decay rate 
from Fig. |10[ taking into account any normalization or 
energy scaling. 

Our analysis uses simple data reduction, minimal back¬ 
ground assumptions, and straightforward analysis proce¬ 
dures. Thus there are many ways our results can be 
improved with further study. 

Firstly, the dominant uncertainty on the data currently 
comes from the energy behavior of the effective area. In 
this work, we use a constant 5% uncertainty as quoted by 
the GBM collaboration, which we then validate against a 
power-law assumption. The limit can be improved if this 
uncertainty can be better quantified or calibrated, e.g., 
using known background lines. 

Secondly, most of the observed counts come from 
cosmic-ray related events. In this work we do not at¬ 
tempt to model such background in detail. In principle, 
this background can be better understood with simula¬ 
tions of cosmic-ray interactions with the satellite that 
take into account the satellite geometry and composi¬ 
tion. Additionally, one can characterize the cosmic-ray 
induced background by satellite positioning, either by 
using high energy observations where the data is back- 
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FIG. 11: Constraints from X-ray missions on sterile neu¬ 
trino dark matter decays, which depends on the mixing angle, 
sin 2 (26), and the mass, m a . In the shown mass range, the best 
previous constraints are set by observations of cosmic X-ray 
background (CXB) from HEAO-1 (50] and the Milky Way 
(MW) halo from INTEGRAL jB(J. The lower bound of the 
parameter space (black-hashed region) is valid for the model 
i/MSM [33]. The flux (spectral) analysis limit derived from 
this work is shown in black-dashed (blue) line. 


ground dominated, or using the Earth’s magnetic field 
information. One can possibly construct a template for 
cosmic-ray induced events, which allows background re¬ 
duction using spatial information. A significant reduc¬ 
tion of the cosmic-ray background can further improve 
the dark matter limit for all energies. 

The next source of backgrounds for dark matter 
searches is those arising from astrophysical origins, which 
dominate at low energies. This background can be mod¬ 
eled using high resolution X-ray sky maps from other 
missions. One can generate an astrophysical template 
for the GBM detectors which can then be used to sub¬ 
tract the astrophysical contribution from the data. 

Importantly, if eventually the data are reduced to a 
regime where statistical uncertainties become important, 
it is important to treat systematic and statistical uncer¬ 
tainty simultaneously mi- 

ill principle, the analysis can be extended to higher 
energies. The GBM-Nal detectors are sensitive up to 
1 MeV and the GBM-BGO detectors extend to 40 MeV, 
which is even higher than INTEGRAL and complemen¬ 
tary to COMPTEL. However, at higher energies, the Nal 
detectors start to observe photons from backward direc¬ 
tions due to either re-scattering or penetrating photons. 
The BGO detectors are designed to be sensitive to both 
front and back directions. A dedicated analysis taking 
into account this detector response for signal modeling is 
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therefore necessary. 

B. Astrophysical Implications 

It is important to highlight that our novel use of GBM 
successfully detects the Galactic astrophysical compo¬ 
nent at 10-20keV energies, as shown in Fig. [ 5 ] and [ 7 ] 
Due to the broad point spread function of the Nal detec¬ 
tor, this seemingly diffuse astrophysical component con¬ 
tains all points source emission near the GC, and possibly 
some diffuse emission. GBM observations may be used 
to impose interesting limits on the total Galactic astro- 
physical emission, which would extend the results from 
INTEGRAL H3H55] down to ~ 10 keV. 

To constrain the astrophysical component, it is nec¬ 
essary to further reduce the detector background (see 
suggestions in the previous section). The analysis pro¬ 
cedures would also need to be modified for a continuum 
spectrum. Such an analysis and a detailed interpretation 
of the astrophysical component are beyond the scope of 
this work. 


C. Conclusion 

We use data obtained by a GBM Nal detector (det-7), 
on board Fermi to set limits on dark matter decaying into 
mono-energetic photons. We first perform a conservative 
flux analysis, based on comparing the total flux normal¬ 
ization of the data and the model. We then perform 
a spectral analysis that assumes the total of non-dark 
matter contributions exhibits a power-law flux spectrum 
within the search window. Our spectral analysis is able 
to improve the limit by about an order of magnitude 
compared to previous searches using CXB with HEAO-1 
data for line energies in the 10-25 keV energy range, or 
20-50 keV sterile neutrino mass. 

Conventionally used to detect and locate transients 
such as gamma-ray bursts, this is the first time that the 


GBM has been used to construct an all-sky map and 
search for dark matter emissions. After performing care¬ 
ful background reduction procedures, we are able to de¬ 
tect the astrophysical component centered on the GC. 

Dark matter searches with GBM benefit from the large 
sky coverage and good energy resolution. Although the 
GBM does not have excellent angular resolution, this is 
not a severe problem for decaying dark matter for which 
the expected Milky Way signal is appreciably extended. 
A unique advantage of the GBM-Nal detectors is that 
they are sensitive to an energy range that is too high for 
X-ray telescopes such as Chandra , Suzaku , and XMM- 
Newton, but too low for INTEGRAL , therefore filling an 
energy gap that was last probed by HEAO-1 in the late 
1970s. 

The current search sensitivity is dominated by system¬ 
atic uncertainties in the effective area. A better under¬ 
standing of GBM detector response through simulation 
or calibration, as well as a better understanding of detec¬ 
tor and astrophysical backgrounds, can further improve 
the data quality and resulting limits substantially. 
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